The Fermi accelerator in atom optics 
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We study the classical and quantum dynamics of a Fermi accelerator realized by an atom bouncing 
off a modulated atomic mirror. We find that in a window of the modulation amplitude dynamical 
localization occurs in both position and momentum. A recent experiment [A. Steane, P. Szriftgiser, 
P. Desbiolles, and J. Dalibard, Phys. Rev. Lett. 74, 4972 (f995)] shows that this system can be 
implemented experimentally. 
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Almost fifty years ago Enrico Fermi |lj suggested that 
it is collisions with moving magnetic fields that accelerate 
cosmic rays. Since then many models associated with the 
name of Fermi accelerator describing particles colliding 
with moving walls have been investigated theoretically 
The recent experiment on atoms bouncing off 
a modulated atomic mirror constitutes an atom optics 
H version of the Fermi accelerator. In the present pa- 
per, we argue that this system serves as a testing ground 
for classical and quantum chaos || . We show that there 
exists an experimentally accessible parameter regime in 
which such experiments can bring out many features of 
the Fermi accelerator. 

Three properties make this model rather special |T^] 
in the field of chaos: (i) It is one of the very few bound 
systems with no continuum pd] , p^ |, (ii) the phenomenon 
of dynamical localization occurs in the position and mo- 
mentum variables [fl3| , (iii) dynamical localization arises 
only in a window of the modulation depth; the onset of 
classical chaos sets the lower boundary, whereas the ef- 
fective dimensionless Planck constant determines the 
upper boundary. Therefore, we can tune the system con- 
tinuously from a regime of no chaos, through one with 
classical chaos but dynamical localization, to one with 
classical chaos and no dynamical localization. The ex- 
periments performed so far have focused on the latter 
regime. However, we show that they can easily be ex- 
tended into the localization window. 

We now consider a cloud of laser-cooled atoms stored 
in a magneto-optical trap. When we switch off the trap, 
the atoms move along the 5-direction [ pd[ under the influ- 
ence of gravity and bounce off an atomic mirror ]l5| ] . The 
latter results from a laser field incident on a glass prism 
under an angle of total internal reflection. This creates 
an evanescent wave whose intensity I(z) = Iq exp(— 2kz) 
decays over a distance fc _1 outside of the prism. A 
sinusoidal modulation of the intensity with ampli- 
tude e and frequency u> changes the intensity I(z) to 
I(z,t) = Io exp(— 2kz)(l + esmujt). In our calculations 
we have used a slightly different form 



I(z, t) — Iq exp(— 2kz + e sinwi) 



(1) 



mirror and is our exponential model of the Fermi acceler- 
ator. When e is not too large, the results do not depend 
significantly on the form of the modulation. 

In order to avoid problems associated with sponta- 
neous emission we consider a large detuning between the 
laser light field and the atomic transition frequency. This 
ensures that the atom rarely leaves the ground state. The 
dynamics of the center-of-mass motion of the atom in the 
ground state follows then from the Hamiltonian 



P 2 

H = h mgz 

2m 



eff e ~2kz+e sin Lut 



(2) 



Here p is the momentum of the atom of mass m along 
the z-axis, and g denotes the gravitational acceleration. 

We introduce the dimensionless position and momen- 
tum coordinates z = zuj 2 /g and p = puj/(mg) and 
time t = cot together with the dimensionless intensity 
Vo = fiuj 2 n, e f f I (4mg 2 ), steepness k = 2kg/uj 2 and the 
modulation depth A = u> 2 e/(2kg) of the evanescent wave. 
In these variables, the classical dynamics follows from the 
Hamilton equations of motion, 



z = p, 

p = —1 + kVq exp [— k(z — Asini)] 



(3) 



The corresponding Schrodinger equation for the atoms 
in the ground state reads 



9-0 



p 

— + z + Vq cxp [— k(z — A sin £)] 



V>, (4) 



of modulation which corresponds to an oscillation of the 



where & = hto 3 / (mg 2 ) denotes the dimensionless Planck 
constant Jig ], 

In this paper, we use Eqs. (||) and ([I]) to determine 
the classical and quantum mechanical position and mo- 
mentum distributions of cold atoms bouncing under the 
influence of gravity off an oscillating mirror. In all our 
calculations we start at t = from an ensemble of atoms 
with an average momentum zero and an average position 
z = 20 above the mirror. We find dynamical localization 
in the quantum case. However, localization occurs only 
over a certain range of the modulation depth A. In order 
to understand this result and to find the lower and upper 
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boundaries A; and A u of this regime, we approximate the 
exponential potential by a hard wall, that is an infinitely 
steep and infinitely high wall and borrow some results 
P,|6|JT7|] for such a time dependent triangular potential 
well. 

We can approximate the dynamics of this system by a 
map connecting two consecutive bounces of the atom: 
The momentum p of the atom and the phase 9 of the 
wall before a bounce determine the momentum p and 
phase 9 just before the next bounce through the stan- 
dard Chirikov- Taylor map 

p = p + K cos 9 

9 = 9+p mod(27r) (5) 

with K = 4A. 

Chirikov |1| has shown that when the chaos parameter 
K becomes larger than the critical value K cr = 0.9716..., 
the classical system undergoes a global diffusion. This 
implies for our driven triangular potential well that diffu- 
sive dynamics sets in for A > A; = K cr /4 w 0.24. Below 
this value we have isolated resonances, whereas above, 
the resonances overlap and we have islands embedded in 
a stochastic sea. 

In the corresponding quantum mechanical system, the 
quasi-energy spectrum of the Floquet operator changes 
from a point spectrum to an almost continuum || when 
A > \ u = Vfc /2. Here quantum diffusion destroys local- 
ization. The conditions of classical and quantum diffu- 
sion, together, define the window 

0.24 < A < -y- (6) 

in the modulation strength, where we can find dynamical 
localization. 

This result obtained for the triangular well determines 
an approximate range of A in which we can observe dy- 
namical localization for the exponential well. Since the 
lower boundary A/ is set by classical dynamics we can find 
A; by evaluating Lyapunov exponents. For a modulation 
amplitude A < 0.24, the Lyapunov exponent converges to 
zero in a vast range of initial conditions except in small 
regions near separatrices as shown in Fig. 0(a). However, 
for larger modulations diffusion occurs and the Lyapunov 
exponent becomes positive in a vast range of initial con- 
ditions as shown in Fig. |l|(b). 

Above the lower boundary, the classical system un- 
dergoes diffusion in both position and momentum space. 
This diffusion manifests itself in a linear growth of the 
square Ap 2 of the width Ap = ((p 2 ) — (p) 2 ) 1 / 2 of the 
classical momentum distribution and of the width Az = 
((z 2 ) — [z) 2 ) 1 / 2 of the classical position distribution with 
time as shown in Fig. |^. The linear growth of Ap 2 with 
time, Ap 2 ~ Dt, also follows from the Chirikov- Taylor 
map H of Eq. ([)]). However, the value of the diffu- 
sion constant obtained for our choice of parameters is 



much smaller than the value D = K 2 /2 predicted by the 
Chirikov- Taylor map. 

We note that the numerical result of Fig. || also sug- 
gests a linear time dependence of the width Az of the po- 
sition distribution. In addition we see that the values of 
Ap 2 and Az are approximately equal. These facts can be 
explained by assuming that the distribution of positions 
and momenta is governed by the classical Boltzmann dis- 
tribution 

P cl (z,p) = (2n)- 1 / 2 rr 3/2 exp[-(p 2 /2 + V(z))/r ] ] , (7) 
where 

V(z) = z + V exp(-Kz). (8) 

Here, the quantity r\ depends on time and plays the role 
of an effective temperature. This conjecture is supported 
by the fact that the Boltzmann statistics holds in case of 
diffusive dynamics ||l9|| . 

The calculation of the Az and Ap 2 can be performed 
analytically in the simplest case of the triangular poten- 
tial well and it gives the equality 

Az = r] = Ap 2 . (9) 

We have checked that for an exponential barrier this 
equality is also approximately true. 

We also show in Fig. [}] the average momentum p and 
average position z. Classically, the average momentum 
oscillates around zero, which corresponds to our initial 
average momentum, in agreement with the result ob- 
tained from Chirikov map. On the other hand, the av- 
erage position displays a linear rise with time, in accor- 
dance with the relation calculated with the Boltzmann 
distribution. 

The classical position and momentum distributions 
shown in Fig. ^ follow then the exponential barometric 
formula P c i(z) — 77 exp(— z/rj) and the Gaussian distri- 
bution Pciip) = (2-kt))~ x I 2 exp[— p 2 /(2rj)] predicted from 
Eq. (0). 

With the help of these distributions we can easily es- 
tablish the relation Az = 77 = Ap 2 . Since we have the 
diffusion law Ap 2 ~ Dt, Eq. (||) predicts that Az ~ Dt 
in agreement with the numerical results of Fig. |2| 

In contrast to this classical diffusion the corresponding 
quantum mechanical quantities saturate after an initial 
rise that is of classical nature. This difference between 
classical and quantum dynamics manifests itself after the 
quantum break time t*. We estimate [B| this time as 

— 2 . 

t* ~ D"k and the corresponding saturation value of 
Ap 2 follows from Ap 2 ~ Dt* ~ D 2 1t~ 2 . 

We also note characteristic oscillations in the quantum 
mechanical widths Ap 2 and Az. These oscillations, ab- 
sent in the classical curves, are a generalization to driven 
systems plj of the revival phenomena |^0| . 
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The quantum distributions in position and momentum 
shown in Fig. [| are completely different from their clas- 
sical counterparts. Indeed the quantum mechanical mo- 
mentum distribution is exponential rather than Gaus- 
sian. Moreover, the quantum mechanical position distri- 
bution contains two exponentials: The steep one corre- 
sponds to dynamical localization in position whereas the 
flat one has the same steepness as the classical baromet- 
ric formula. This separation of the quantum distribution 
into a quantum and a classical part also occurs in other 
bound problems ]l2] ). 

A qualitative comparison between the quantum trian- 
gular well and the kicked rotator model || yields the 
localization length as I ~ Dk ~ 8A 2 /£ . In our nu- 
merical investigations we clearly see this qualitative de- 
pendence on and A. 

The initial minimum uncertainty wavepacket delo- 
calizes when the modulation depth exceeds the upper 
boundary X u of the localization window. In Fig. we 
show the widths of the classical (thick line) and the quan- 
tum mechanical (thin line) momentum distribution as a 
function of the modulation amplitude A. We find ex- 
ponential localization within the window, as shown in 
the inset (a). However, above the upper boundary A u , 
quantum diffusion sets in resulting in the Gaussian dis- 
tribution shown in the inset (b). 

We now make contact with the recent experiment 
JjJ on the phase modulation of de Broglie waves. In 
this experiment cesium atoms bounce off a mirror with 
fc _1 = 0.19 /im. The mirror is modulated with an am- 
plitude e = 0.82 and various frequencies of the order 
u « 2vr x 900 kHz. With the mass m = 2.21 x 10~ 25 kg of 
Cs atoms, the gravitational acceleration g — 9.81 m/s 2 , 
and fi = 6.673 x 10~ 34 /27r Js, we find the dimensionless 
Planck constant % w 9 x 10 8 and the modulation depth 
A w 2.5 x 10 5 . We emphasize that this value is larger 
than X u = Vk/2 ~ 1.5 x 10 4 . Therefore this experiment 
lies outside of the localization window. 

However, a modulation frequency of w = 2ir x 
1.477 kHz, the decay length fc _1 = 0.455 /Ltm and the 
effective Rabi frequency f2 e // = 2n x 88.8 kHz lead to 
£ = 4, k = 0.5 and Vq = 60 which are the values used 
throughout our paper. We have chosen this value of Vq to 
guarantee that the atoms will not hit the surface of the 
mirror which is situated at z — 0. Of course some atoms 
have enough energy to break through the barrier but 
their number is negligible. The lower and upper bound- 
ary A; = 0.24 and X u — 1 then translate into an intensity 
modulation of e; = 0.12 and e u = 0.5, respectively. 

With our parameters we have observed dynamical lo- 
calization within 100 bounces. The quantum break time 
t* ~ 75ms corresponds to 50 bounces. Since Ref. |22| re- 
ports more than hundred bounces, this effect should be 
observable. 

Note that the initial condition in M was zq = 3.3 mm 



which corresponds to zq — 2.9 x 10 4 in dimensionless co- 
ordinates. In contrast, in our numerical calculations we 
have chosen zq — 20 which corresponds to z = 2.27 /im. 
However, a recent experiment p2| shows that even this 
is possible. We are aware of the fact that the practical 
implementation of our proposed experiment may not be 
trivial. However, in spite of the experimental challenge — 
due to the external modulation of the mirror — we are 
confident that an experiment along the above lines can 
be performed. 

We conclude by summarizing our main results. An 
atom bouncing off a modulated mirror under the influ- 
ence of gravity exhibits dynamical localization in posi- 
tion and momentum. However, this effect only occurs in 
a window of modulation. Our investigations show that 
bound systems may exhibit a rich dynamical behavior 
both in classical and quantum domain, which is substan- 
tially different from the standard kicked rotator model. 
These conclusions motivate further theoretical and ex- 
perimental studies of bound systems from the view point 
of chaos. Since our system in the absence of the driv- 
ing force does not contain any continuum of states it is 
cleaner than the microwave driven hydrogen atom. Cur- 
rently available experimental technology allows us to ob- 
serve this phenomenon. 

Our analysis of the atomic Fermi accelerator is based 
on a laser mirror for the atom and therefore makes use 
of the interaction of the atom with an evanescent laser 
field. However, it is interesting to note that also mag- 
netic mirrors p3| for atoms exist and have produced 
many bounces. A modulation of such a magnetic mir- 
ror could offer another possible realization of the atomic 
Fermi accelerator. 
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FIG. 1. Lyapunov exponent C of an atomic Fermi accelerator below (a) and above (b) the lower boundary A; = 0.24 of 
the localization window. We use three initial conditions (zo,Po) = {(20, 0); (20, —2); (40, —2)} represented by thick dots in the 
Poincare sections. In (a) we have chosen the modulation depth A = 0.2 and all initial conditions lie inside isolated resonances. 
Consequently for all three initial conditions the Lyapunov exponent approaches zero. In (b) we have A = 0.5 > A;. Here, the 
phase space point (40,-2) still sits in an island whereas the points (20,0) and (20,-2) lie in the stochastic sea. As a result the 
Lyapunov exponent for the first initial condition converges to zero whereas for the other two it is positive. Here and in all other 
figures we have chosen for the height of the potential Vo = 60 and for its steepness k = 0.5. 
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FIG. 2. Comparison between quantum (thin lines) and classical (thick lines) values of the (a) square of the width Ap of 
the momentum distribution and (b) width of the position distribution Az. In the quantum case we integrate the Schrodinger 
equation, Eq. (Q) , subjected to the initial condition of a Gaussian minimum uncertainty wavepacket located at the phase space 
point (20,0) with a width Az — 2 in position and the corresponding width Ap = ~k/(2Az) in momentum. In the classical case 
we propagate an ensemble of particles distributed according to the same Gaussian distribution using the classical Hamilton 
equations, Eq. (^|). The number of particles in the classical simulation is 2000. The height of the exponential potential is 
Vo = 60, its steepness k — 0.5 and the modulation strength is A = 0.5. In the quantum mechanical case the effective Planck's 
constant is taken as % — 4. 
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FIG. 3. Comparison between the classical (thick lines) and quantum (thin lines) average momentum (a) and average position 
(b) as functions of time. For the sake of presentation, we show the envelope of the corresponding functions. We note the 
modulation of the quantum mechanical envelopes which is absent in the classical curves. This is a manifestation of revivals |2(]] 
in a driven quantum system [Ell. The calculations were performed for the same set of parameters as in Fig. 0. 
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FIG. 4. Comparison between the quantum mechanical (thin lines) and classical (thick lines) momentum (a) and position (b) 
distributions on a logarithmic scale for an atomic Fermi accelerator. The momentum distribution in the quantum mechanical 
case exhibits an exponential localization, whereas the corresponding classical distribution is Gaussian. The classical as well as 
the quantum mechanical position distributions are both of exponential form. However in the classical case this form results 
from the linear potential in the Boltzmann distribution. The peak around z — 20 is due to the fact that a considerable part of 
our initial ensemble lies inside a stable island. In the quantum case we find two exponentials: the flat one is a remnant of the 
classical Boltzmann distribution and the steeper one represents dynamical localization. The dashed lines indicate linear and 
quadratic fits which correspond to exponential and Gaussian distributions, respectively. Here we have chosen A = 0.8 and the 
number of particles in the classical simulation is 10000. The integration time is t — 2650 and all the other parameters are the 
same as in Fig 0. 
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FIG. 5. Localization window denned by the onset of classical and quantum diffusion at A; = 0.24 and \ u = sfk/2 = 1, 
respectively. Squares of the width of the classical (thick line) and the quantum mechanical (thin line) momentum distributions 
in their dependence on A. The two curves start to separate at A = A;. A transition from a localized to a delocalized quantum 
mechanical momentum distribution occurs at A„. For A = 0.8 which for ~k = 4 lies well within the window we find exponential 
localization [inset (a)]. In contrast for A = 1.2 which lies outside of the window we find a broad Gaussian distribution indicating 
derealization [inset (b)]. All the other parameters are as in Fig. | and t = 3200. 
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